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Abstract 

We calculate, in a class of Gauge invariant functionals, by variational meth- 
ods, the difference of vacuum energy between two different backgrounds: 
Schwarzschild and Flat Space. We perform this evaluation in an Hamilto- 
nian formulation of Quantum Gravity by standard "3 + 1" decomposition. 
After the decomposition the scalar curvature is expanded to second order 
with respect to the Schwarzschild metric. We evaluate this energy difference 
in momentum space, in the lowest possible state (regardless of any negative 
mode). We find a singular behaviour in the UV-limit, due to the presence of 
the horizon when r = 2m. When r > 2m this singular behaviour disappears, 
which is in agreement with various other models presented in the literature. 



I. INTRODUCTION 



The problem of computing quantum corrections to a classical energy in a complicated 
theory such as Einstein gravity, can be approached by performing an analysis of the thermo- 
dynamical quantities that characterize the system under consideration. This analysis can 
be done by means of the computation of the free energy of the system at a given volume 



1 



and temperature. Defining the Euclidean action as 

1 M = -J^r L d * x ^ R ^ -^rL d " x ^ K ^ C 1 ) 

where i? is the Ricci scalar of the metric g^ v and K\ is the trace of the second fundamental 
form, one can compute quantum corrections to the Euclidean action, of a fixed background 
geometry with the appropriate boundary conditions. Since we wish to study quantum fluc- 
tuations with respect to the Schwarzschild geometry, there are two types of boundary 
conditions that are related to the background under consideration: 

a) Asymptotically Euclidean (AE), 

b) Asymptotically Flat (AF). 

An AE background metric is one in which the metric approaches the flat metric on i? 4 
outside some compact set. The boundary at infinity is topologically S 3 . An AF background 
metric is one in which the metric approaches the flat metric R 3 x S 1 outside some compact 
set. The boundary of infinity is topologically S 2 x S 1 . Anyway, we could consider a different 
point of view based on the Hamiltonian approach. In this framework one is able to deal with 
three dimensional fields configurations separated out by the time variable. The advantage 
of the Hamiltonian framework is that one can manage from the beginning with energy fields 
configurations which give directly the quantum corrections to the classical term. To do this, 
the first step is the separation of the three dimensional space from the time by means of the 
ADM. variables Iffl. In terms of these variables the line element becomes 



ds 2 = g^ (x) dx^dx v = —N 2 (dx°) + g iS (nUx + dx 1 ) (N j dx° + dx j ) 

-N 2 + NiN*) (dx°Y + 2Njdx°dx j + g ij dx i dx j . (2) 



iV is called the lapse function, while iVj the shift function. The associated matrix represen- 
tation of g^u is 



I -N 2 + N i N i Nj^ 
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(3) 



with the inverse given by 



9 : 



1 

N 2 



N J 
N 2 



\ n 2 y % 3 



TV, TV 1 



(4) 



N 2 j 



Roman indices will be raised and lowered by the induced metric on the three surface x°. 
In terms of the ADM. variables, the initial action can be written as a sum of a "kinetic" 
and a "potential 1 term 



^ Jdx°N J dx 3 ■\<7]{(K IJ K« - K 2 ) + Mr} , 



(5) 



where = ^ + Nj\i — fi^o) is called the second fundamental form and "| " means 
covariant differentiation with respect to the 3D gravitational background, K is the trace of 
the second fundamental form, 3 R is the scalar curvature in 3D, ^ y/g is the invariant of the 
metric in 3D. In this form the time derivative is isolated and it is possible the computation 
of the conjugate momentum to gij, that is 

SI 



TV 



Kl 



By a Legendre transformation we calculate the Hamiltonian 

1 



(6) 



H 



d 6 x 



d 3 xlN 



16nG 



KijK ij — K 2 ) + (3) R N (3) ^ 



16nG 



TTijir 



7T2 



(3) 



+ Ni(2TT 



(7) 



Wy/g V" 1 " 2 ) v ^167rG 

The first term of (0) has a quadratic structure in the momenta, suggesting, as a first approx- 
imation that we could compute quantum corrections to the energy expanding 3 R in terms 
of the quantum field fluctuations with respect to a given background, e.g. the Schwarzschild 
background. Since we wish to understand the pure gravitational vacuum, we neglect the 
matter fields and since in our approach only the spatial part of the background comes into 
play our background is of the wormhole type [jnj. Although the energy computation at 
quantum level is very unclear because of the constraint coming from the lapse function we 
will adopt the Hamiltonian approach the same and to this purpose a simple framework will 



be shown in section [H[ The rest of the paper is structured as follows, in section [III] we 



analyze the orthogonal decomposition of the Hamiltonian both in tangent and co-tangent 
space, in section [TV] we define the gaussian wave functional for gravity in analogy with non- 
abelian gauge theories, in section [V] we give some of the basic rules to perform the functional 
integration and we define the Hamiltonian approximated up to second order, in section [VI|, 
we analyze the spin-2 operator or the operator acting on transverse traceless tensors, only 



for positive values of E . We summarize and conclude in section [VTi 



II. THE HAMILTONIAN ON THE SLICE 



After the introduction of AVA4 

variables, we recall that the Hamiltonian is: 



H = f d 3 x(NH+N i H i ) (8) 



where 



H=G ijkl Tr ij Tr kl (J^j - (^fj (3) i? (Super Hamiltonian) (9) 

and 

TC = -27rg (Super Momentum). (10) 

In ([!(]) the derivative is covariant with respect to the 3D background field, ti is the usual 
Planck mass, and Gijki is the Wheeler-DeWitt (WDW) metric. If we look at N and Ni as 
fundamental objects describing the correct variables, by variational principles we obtain the 
usual constraint equations, that is 

U =0,W =0 Classical 

(11) 

Hty = 0, W * = Quantum 

The usual interpretation of these equations is that they represent constraints on the initial 
value problem or in other words they represent gauge invariance with respect to time and 



gauge transformations. Nevertheless we have a chance to define and computing energy if 
we restrict on a given hypersurface, fixing the lapse function to a constant. Such a gauge 
choice is the most appropriate for wormhole configurations of the background geometry |I[ 
and in particular for the Schwarzschild wormhole. By rescaling time intervals, we obtain^ 

N=l. (12) 

Actually, this choice is compatible with the suspension constraint that one can adopt in 
quantum cosmology to obtain a Schrddinger-like equation, provided at the end of the cal- 
culation one assures that the gauge invariance is restored [|TJ ||. Then the Hamiltonian in 
the time-like gauge is 



H = / d 3 xH = / d 3 x 



(13) 



Since (|13|) is valid on a " fixed 1 hypersurface, to recover the original equation, i.e. (|8[), the 
correct procedure to perform will be a summation over all possible lapses; this means that 
the constraint (|TT|) (at classical or at quantum level) will be restored after this summation. 

III. ULTRALOCAL METRICS AS A TOOL FOR DECOMPOSING TENSOR 

FIELDS 

Instead of performing calculations in the usual WDW metric we will use a one-parameter 
family of supermetrics to disentangle gauge modes from physical deformations. For this 
reason we require an orthogonal decomposition for both 7^ ■ and , that is we need a metric 
on the space of deformations, i.e. a quadratic form on the tangent space at h. The condition 
of ultralocality, where Gijki locally depends on but not on its derivatives, could be taken 
as a good condition for the functional measure, explicitly: 



1 A different treatment, but close to our present approach, is based on the separation between 
dynamical variables and embedding parameters and can be found in Ref. [Hj 
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(h,k) :-- 



M 



J-gGf% (x) k kl (x) d 3 x, 



(14) 



where 



Qijkl = ( g ik g jl + g il g jk _ ^OLQ^Q 



Id \ 



(15) 



The WDW metric, introduced in (|9]), is just ( 15 ) with a = 1. The "inverse" metric is defined 
on co-tangent space and it assumes the form 



(p,q):= I ^G^(x)q H (x)d 3 x, 



M 



(16) 



where 

Gijki = (9ik9ji + 9ii9 jk ~ W9ij9ki)- ( 17 ) 

with a + (3 = 3a/3, so that 

GT m Gn mk i = \{^l + 5\5{). (18) 

These are non-degenerate bilinear forms for a ^ for a = j- the metric is not invertible 
and becomes a projector onto the tracefree subspace, while is positive definite for a < gr- 
and of mixed signature for a > j- with infinitely many plus as well as minus signs. 

We have now the desired decomposition on the tangent space of 3-metric deformations 

h v = \h9i 3 + (Li) i3 + h^ (19) 

or, in matrix form, 

h = ^hg + (RangeL) + (KerL ] ) , (20) 

where the operator L maps ^ into symmetric tracefree tensors, according to @ ||, 
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(LO tJ = V£ + V£ - \g l3 (V • • (21) 

Consequently, the inversion of the metric ([H]) (that is fllTD) guarantees us the same decom- 
position also in phase space (co-tangent space). 

IV. THE GAUSSIAN WAVE FUNCTIONAL 

There are some reasons to introduce a gaussian wave functional for the description of 
the vacuum state in gravity. Starting from the analogy between nonabelian gauge theories 
and gravity we illustrate how gaussian wave functional work in the former case. We define 
1 

v& [At (If)] = TVexp |-i J d 3 xd 3 y5A« (if) (if, y>) M{ (V)} (22) 
where A/" is a normalization factor and where 

^(lf) = ^(^)-Af(^). (23) 
In equation fl23|) , A" ("of) is a background field which can be treated as a variational param- 



eter together to the function Gf- (if, if) in fl22|) . 



From the definition in (^) one finds the expectation values 

(^|4 a (x>) = (x>) , 

(*| A ? (^) 4 (V) |*) = 3? 0?) ^ (IT) + g# (if, V) , 

(*|£?C^)|*) = 0, (24) 
(if) Ej (iT) I*) = (if, if) , 

(If) |*) = W(^) + K-fc/^Gg- (If, If) , 
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where 



Ef 



—i 



8 At (it) 



(25) 



and 



Bt (*) = U jk [V,AI (it) + \]^A\ (1?) A) (f)) 



(26) 



eijk is the usual anti-symmetric tensor and / are the structure constants of the gauge 
group, for ex. SU (N). After having experienced how the apparatus works on nonabelian 
gauge theories, we define a "Vacuum Trial State" for gravity, and for this purpose we recall 
the orthogonal decomposition fll8|) to look at the essential structure of the inner product 
between three-geometries 



(h, h) := / JgGf% (x) h u (x) d*x 



M 



M 



- - a) h* + (L£f {L$ u + h^h 



(27) 



Previous formula leads us towards the definition of the "possible" trial wave functional for 
the gravitational ground state 



hu 0^) = A^exp < — 



4/H 



hK- i h\ ± +(m)K- i m)\ +(hK- i h) 



Trace 



I x,y 



, (28) 



or in other terms 



h 3 c*; 



K 3 (it) 



(29) 



In ( p8| ) and in ( p9|) , is the tracefree-transverse part of the 3D quantum field, i s 
the longitudinal part and finally h is the trace part of the same field. The dependence of 
the functional by a will not be discussed in this paper. 

In (|28"D , (-, •) denotes space integration and K~ l is the inverse propagator. The main 
reason for a similar ll Ansat^ comes not only from (|27| ) but even to the observation that the 
momenta quadratic part of the Hamiltonian decouples in the same way. Even if we had to 
give up to ( ^q) from the beginning, making a more general u Ansat£ J about the vacuum wave 



S 



functional (and for more general we mean eqn. (|22|) ) one would discover that the kinetic 
part decouples in these three terms. For completeness, we give the analogous expectation 
values for TT tensors. The other components satisfy the same rules 



^l^(^) m = g. 



x 



*\9i3 (#) 9m (it) I*) = ffti C#) 9m (It) + Ktm 0^ t) 



ijk 



*|7T±(2>)|*) = 0, 



(30) 



where tt, 



'.J 



' Sh-h. (x) 



(*) <i (t) i*) = (t, t) , 



is the representation for the TT momentum. 



V. ENERGY DENSITY CALCULATION IN SCHRODINGER 

REPRESENTATION 



To calculate the energy density associated to the trial functional, we need to know the 



action of some basic operators on \1/ 



. The action of the operator on 



hi Alt) is realized by 



h ii (x)\*) = h ii (x)V \A-at) . 



The action of the operator 7r-- on |^), in general, is 



(x) !*> = -« 



5hij (x) 



The inner product is defined by the functional integration: 



(*l|*2>= / [Vhyix)] *!{/>y}M*«}. 



and the energy eigenstates satisfy the Schrodinger equation: 
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Sh^ (x) 



, hij (x) * [hij (If)} = EV {h^ (If)} , 



(31) 



(32) 



(33) 



(34) 



where H j— i -^ 5 ^ , hij is the Hamiltonian density. Instead of solving fl34|), which is of 
course impossible, we can formulate the same problem by means of a variational principle. 
We demand that 

/,T, I ZJ I ,T,\ f \T>n±(rW r rflrXb* in^Xn^l 

(35) 



<* 




*) _ / 


V9ij (4 


/d 3 xvl/*{^}^vl/{^} 


<* 
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be stationary against arbitrary variations of \I/ j/iy ("a?)}. The form of | H | 1 J r ) can be 
computed as follows. We define normalized mean values by a straightforward modification 
of (BUI), i.e. 



j[Vg^(x)]fd s xg^(x) {,/,• (,•)} 



f[Vg±(x)} \*{g±} 



(36) 



3« 0) 0iu (a:) + V) 







Jd 3 xg^(x)gi l (y) | * {g± (x)} | 2 




^ !•'•) 


I 2 



It follows that 



(x)] (x) - ^ (x)) | * {gjj (*)} | 2 = 
by translation invariance of the measure 



T>H± (x 



hj ] (x)\^{g^(x)+g^(x)}\ 2 =0 
hj- J (x)\^{hj j (x)}\ 2 =0, 



and (|37|) becomes 



/ [Ztti «•'•) / (*) 4 (v) I * [Hi (*)} 



Vhf 3 (x) 



(37) 



(38) 



(39) 



Rather than applying the variational principle arbitrarily, the gaussian Ansatz is made, 
according to which in the beginning of this calculus one has to replace previous general 
formulas with 
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/ij,- (if) = AAexp 



77) K- 1 



+ 



.xvy 



(40) 



With this choice and with formulas (|38|, |39|), the one loop-like Hamiltonian can be written 

as 



A1 2 p Jm 



d'x^gGf 1 [K-^ (x, x) ijkl + ( A 2 )J K± (x, x 



iakl 



(41) 



where the first term in square brackets comes from the kinetic part and the second comes 
from the expansion of the 3 R up to second order in such a way to obtain a quantum harmonic 
oscillator equation type, the Green function K 1 - (x,x) iakl can be represented as 

hi (x) (y) 



K± ^ X ^'-^ 2X N (p) ' 



(42) 



where hj~ a (x) are the eigenf unctions of A£ ■ and (p) are infinite variational parameters. In 
formula (|4~l|) we have written the Spin-2 contribution to the energy density alone; expressions 
like (pi) exist for Spin-1 and Spin-0 terms of 7i. 



VI. THE SPECTRUM OF THE SPIN-2 OPERATOR AND THE EVALUATION 

OF THE ENERGY DENSITY 



The Spin-2 operator is defined by: 



A 9 := -A + 2Ric 



(43) 



or in components, 



:A 9 r : = -a# + 2w 



(44) 



where A is the curved Laplacian (Laplace-Beltrami operator) on a Schwarzschild background 
and 

Rj is the mixed Ricci tensor whose components are: 
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f — 2m mm] . . 

^ = (45) 



where 2m = 2MG. This operator is similar to the Lichnerowicz operator provided that we 
substitute the Riemann tensor with the Ricci tensor. In (f43|) or Ricci tensor acts as 



a potential on the space of TT tensors; for this reason we are led to study the following 
eigenvalue equation 

(-A5] + 2R«) h\ = E 2 h) (46) 

where E 2 is the eigenvalue of the corresponding equation. In doing so, we follow Regge and 
Wheeler in analyzing the equation into modes of definite frequency, angular momentum and 
parity. In this paper we are interested to positive E 2 and low lying states with L = M = 0, 
where L is the quantum number corresponding to the square of angular momentum and M 
is the quantum number corresponding to the projection of the angular momentum on the 
z-axis. For L = 0, Regge- Wheeler decomposition shows that there are no odd-parity 
perturbations at all, therefore: 

2m\~ 1 



h*] en = diag 



H(r)[l- — ) , r 2 K (r) , r 2 sin 2 dK (r) 



^oo 4>) ■ (47) 



The representation fl4"7|) is very useful, because of the decoupling of the components, in fact 

-AH(r)-^H(r) = E 2 H(r) 

—AK (r) + ^rK (r) = E 2 K (r) (48) 

—AK (r) + (r) = E 2 K (r) 

The Laplacian in this particular geometry can be written as 

/ 2m\ d 2 /2r-3m\ d , ln . 

A= 1-— — -5— h-. (49) 



r / dr 2 \ r 2 J dr 

Defining reduced fields, such as: 



ff(r ) = *M iJf(r ) = *M, (so) 
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and changing variables to 



the system 



x = 2m 
becomes 



r r I \ ( I I 1 



-£,h(x)-V(x)h(x) = E 2 h(x) 



(51) 



■£tk{x) + V{x)k(x) = E 2 k(x) 



(52) 



' dx 2 



where 



k (x) + V (x) k (x) = E 2 k(x) 



We note that the new variable is such that 



x ~ r r 



oo V fx) 



a; ~ r — > r V (x) 



3m 



const, 



(53) 



(54) 



where r is the wormhole radius, satisfying the condition r > 2m, strictly. The solution 
of (|52"D , in both cases (flat and curved one) is a Bessel function and precisely the spherical 
Bessel function of the first kind for the L = value of the angular momentum 

2 sin (pr) / 2 sin (pr) 



j (jpr) = p\ 



7r r 



(55) 



7r 

The corresponding Green function for this problem will be 

K( ) _ h(px)j (py) l 
l ' yJ " 2A 4tt 

Substituting (|56|) in fl4l"l) one gets (after normalization in spin space and after a rescaling of 
the fields in such a way to absorb I 2 ) 



(56) 



y 2 roo 



\{P) 



E 2 (p, m) 
\(p) 



(57) 
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where 



E i, 2 (P> m ) = P 2 =F 



3m 



(58) 



Aj (p) are variational parameters corresponding to the eigenvalues for a (graviton) spin-2 
particle in an external field and V is the volume of the system. 

By minimizing ( |57|) with respect to Aj (p) one obtains Aj (p) = [E 2 (p, m)\ 2 and 



E(m,X) = —J2 d P 2^/E? (p, m) with p 2 > 



3m 



(59) 



The total energy in the presence of the background is 

V 1 



Elm) 



2n 2 



- ( dpp 2 (\jp 2 — c 2 + \Jp 2 + c 2 ) where c 2 = — 3- 
2 Jo \ / Tq 



(60) 



For flat space the calculation is essentially the same with the exception of c 2 = 0. Therefore 
the equivalent of ( |60|) in flat space is 

V 1 r°° f 1 — \ 

11/ ^ ^ (61) 



V 1 f°° ( r~ 

£(0) =2^u dpp 



Now, we are in position to perform the energy difference between d60|)and fl6T|). AE (m) up 
to second order in perturbations is 

V 1 foo 



AE lm) 



2n 2 2 



dpp 



p 2 — c 2 + Jp 2 + c 2 — 2Jp 2 



(62) 



We want to evaluate the UV behaviour of (p2|), therefore 
A£(m) = ^f/ °°dpp 3 



1 + ( £ 
1 p 



becomes for p 2 >> c 2 



V c 4 /- 00 dp 



2?H 8 Jo p 



(63) 
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Introducing a cut-off one gets for the UV limitQ 



dp fUx In(^) (64) 



io p Jo x \c 
and A 2? (m) for high momenta can be estimated by the following expression 

At this point we can compute the total energy, namely the classical contribution plus the 
quantum correction up to second order. Recalling the definition of asymptotic energy for 
an asymptotically flat background, like the Schwarzschild one 

e adm = J™ L ^ ik ' j ~ &3>d ^ 

where is the metric induced on a spacelike hypersurface dj\4 which has a boundary at 
infinity like S 2 , one gets, 



2tt 2 \ rg y 16 V 3m / 2vr 2 \ J 16 V 3MG 



(67) 



One can observe that 

AE (m) -> oo when m -> 0, for r = 2m = 2GM (68) 

and 

AE 1 (m) -> when m -> 0, for r /2m = 2GM. (69) 
Remark We would like to explain the reasons that support the results of formula (|6~5). In 



that formula we introduced a particular value of the radius, which behaves as a regulator 



2 It is known that at one-loop level Gravity is renormalizable only in flat space. In a dimensional 
regularization scheme its contribution to the action is, on shell, proportional to the Euler character 
of the manifold that is nonzero for the Schwarzschild instanton. Although in our approach we are 
working with sections of the original manifold to deal with these divergences one must introduce a 
regulator that indeed appears in the contribution of the energy density. 
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with respect to the horizon approach of the potential. The meaning of this particular value 
is related to the necessity of explaining the dynamical origin of black hole entropy by the 
entanglement entropy mechanism and by the so-called "brick wall modeT ||. Indeed, the 
same mechanism is present when one has to regularize entropy by imposing a kind of cut-off, 
that in coordinate space means r > 2m. 

VII. SUMMARY AND CONCLUSIONS 

The trial wave functional approach, by means of Gaussian configurations, led to possible 
calculations of quantum fluctuations of the gravitational field around some fixed background 
geometry. In particular we have studied a spherically symmetric background and by means of 
Birkhoff 's theorem we can claim that our background is of the Schwarzschild type. Since we 
have performed this analysis without any matter contribution and recalling the definition 
of Ref. |10| our result is valid for a Schwarzschild wormhole. However this calculation 



apparatus is entirely based on the possibility of explicitly breaking the invariance under 
reparametrisations, expressed by the gauge fixing flT2|), leading to the conclusion that the 
final result seems depending on the foliation we choose to work. For this reason to restore 
the invariance under reparametrisations, we need to sum on every lapse functionP]. With the 



gauge choice ([12]) , the problem of defining a correct vacuum energy on every slice is well 
posed and the result shows us an intrinsic energy depending only on the dynamics generated 
by 3-surfaces. 
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3 A detailed version of this procedure will be studied in a future paper [11 
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APPENDIX A: CONVENTIONS 



Here we give the conventions for the metric tensor, connections and the curvature tensor: 



1. Background Metric 



9u 



9 ll = l 



ji> 922 = r 2 , #33 = r 2 sin 2 9 



(Al) 



2. Connection 



r 1 

1 ab 



-( 



1 - ?f)r 



r 2 

1 aft 



r- 1 
r -i o 








— sin cos 9 



- (l - ^ ) rsin 2 # 



r- 1 ^ 



r 3 

1 ab 



cot 9 
r-- 1 cot 9 



3. Riemann tensor, Ricci tensor and the Scalar Curvature in 3D 



(A2) 



R\ 3m = - rj i)m + rj a r^ - rj^r^ Riemann tensor 

Because of the vanishing of the Weyl tensor in 3D, that is C'- m = 0, Riemann tensor 
is completely determined by Ricci tensor 

Rlijm = 9ljRirn ~ 9lmRij ~ 9ijRlm + 9imRu 



R im = R ilm Ricci tensor 



R = g 1 -' Rij Scalar curvature 



17 



APPENDIX B: SCALAR CURVATURE EXPANSION 



In this part we give the necessary tools for the scalar curvature expansion in terms of the 
fluctuations of 3-surfaces around the background geometry. Metric tensor will be separated 
in a classical part (background) plus a quantum part i.e. 

9ij = 9ij + h ij ( B1 ) 



1. Expansion of the determinant 

^fq— = exp Tr In = exp Tr\ In (g^ + h^) = exp Tr\ In g {j (l + |f ) 



expTri [ln^ + ln(l + !£) 



(B2) 



2. The inverse metric expansion 



9 k 



si + K 



gij _ h a + ^ _ + h ik h l k hThP m + O (/> 5 ) (B3) 



3. Connection 



To 0" 1 order r£. = T 4J 



- 2# lSfe,j + - 9ij,l) 
To 1 st order r* = 5* = (h Hlj + h m - /^) 



The higher order corrections to the connection are related by the formula 



r *(n) = _ h k iT Kn-l) where pfc(l) = 5 * 



(B4) 



4. Riemann Tensor 



A- 



p£ p/ , p/ -pa p/ pet 

mij ji,m ' ja mi ma ji 
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It is convenient to divide Riemann tensor into two terms: linear and non-linear 



Lin (R\ jm ) - L\ jm :- T l mi j - T l ji r 



N - Lin (R\ jm ) = N{ 



pi pa pi pa 

ja mi ma ji 



(B5) 



The higher order terms of L\- m are simply 



ri(n) -pl( n ) p'( n ) 

ijm ■ mi\j ji\m 



(B6) 



while higher orders of N l ^ m are 



n-1 



ArZ(n) ._ \-na(j)yl(n-j) _ yl(n-j)-pa(j) 

ly ijm • / j y- mi ja 1 ma ji 



"jjrrt • / j 
3=1 



(B7) 



5. Second order scalar curvature 



Collecting together previous expansion formulas we can write the following expression 
for the 3 R scalar curvature expanded up to second order: 



/ 



<Tx 



--hAh + -h H Ah H - -h i] ViV.ih\ + -hV l V i h H - -h lj R ia h a q + -hR^ 



(B8) 
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